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Abstract 

Based on a previously proposed thermodynamic analysis^ , we study the relative 
stabilities of five Si phases under unaxial compression using ab initio methods. The 
five phases are diamond, /3Sn, simple- hexagonal (sh), simple-cubic, and hep structures. 
The possible phase-transition patterns were investigated by considering the phase tran- 
sitions between any two chosen phases of the five phases. By analyzing the different 
contributions to the relative phase stability, we identified the most important factors 
in reducing the phase-transition pressures at uniaxial compression. We also show that 
it is possible to have phase transitions occur only when the phases are under uniaxial 
compression, in spite of no phase transition when under hydrostatic compression. 

Taking all five phases into consideration, the phase diagram at uniaxial compression 
was constructed for pressures up to 20 GPa. The stable phases were found to be dia- 
mond, /3Sn and sh structures, i.e. the same as those when under hydrostatic condition. 
According to the phase diagram, direct phase transition from the diamond to the sh 
phase is possible if the applied uniaxial pressures, on increasing, satisfy the condition 
of P x > P z . Similarly, the sh-to-/3Sn transition on increasing pressures is also possible 
if the applied uniaxial pressures are varied from the condition of P x > P z , on which the 
phase of sh is stable, to that of P x < P z , on which the /3Sn is stable. 
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1 Introduction 



Most theoretical studies of compressive phase transitions are restricted to hydrostatic com- 
pression. However, nonhydrostatic compression occurs commonly in realistic situations and 
can lead to very different phase-transition pressures. For example, previous experimental 
studies have shownll' HI that nonhydrostatic compression effectively reduces the phase- 
transition pressure of diamond structured Si to the metallic /3Sn phase. The generally used 
enthalpy function H = E + PV in ab initial approach, which has been very successful in 
predicting hydrostatic phase-transition pressure@, is no longer the appropriate quantity 
in determining the relative stability between phases under nonhydrostatic conditions. In a 
previous paper (called paper I hereafter), we have carried out a thermodynamic quantity 
analysis, similar to the use of the enthalpy difference in hydrostatic cases, in order to distin- 
guish the relative stability between phases under uniaxial compression. That was applied to 
the uniaxial phase transitions between the diamond and (3Sn structures of Si and Ge using 
ab initio calculations. It was found that under uniaxial compression the phase-transition 
pressures for both Si and Ge were greatly reduced, i.e. from 11.4 to 3.9 GPa for Si and 9.5 
to 2.5 GPa for Ge. 

In the present work, based on the previous thermodynamic analysis, we study the 
relative stability of five Si phases under uniaxial compression. The five phases are diamond, 
/5Sn, simple hexagonal (sh), simple cubic (sc) and hexagonal closed packed (hep) structures. 
We investigate the possible phase-transition patterns at uniaxial compression by considering 
the phase transitions between any two chosen phases of the five structures. We shall 
show that it is possible to have phases which do not undergo phase transition when under 
hydrostatic compression but do undergo phase transition when under uniaxial compression. 
By analyzing the different contributions to the relativie phase stability, we identify the most 
important factors in reducing the phase-transition pressures at uniaxial compression. We 
demonstrate that the zero-pressure structural parameters of two phases can be used to 
guess the phase-transition pattern of the two phases under uniaxial compression. Finally, 
the phase diagram of the five Si structures under uniaxial compression is presented. 
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The paper is organized as the followings: in the next section the ab initio methods 
used in the studies are specified. The third section describes the hydrostatic results. The 
fourth section presents the ab initially obtained phase-transition patterns at uniaxial com- 
pression and discusses the important factors leading to different phase-transition patterns. 
In the fifth section we present the uniaxial phase diagram of the five Si phases and finally 
the conclusions in the sixth section. 



2 Calculational Methods 

In our studies, the total energy (with respect to the atomic energies) of systems were 
calculated with the ab initio density functional method!!, and the generalized gradient 
approximation (GGA)i for the exchange and correlation functional was used. The calcu- 
lations were performed using the ab-initio total-energy and molecular-dynamics program 
VASP (Vienna Ab-initio Simulation Program) developed at the Institut fur Material Physik 
of the Universitat Wien BIS. Norm-conserving pseudopotentialslMJ were used for the 
electron- ion interactions. All present calculations were carried out in momentum space 
and the Kohn-Sham wavefunctionsEH were expanded in a plane-wave basis. The energy 
cutoffs for the plane-wave basis was 400eV. The integration over the Brillouin zone was 
approximated by k-point sampling using the Monkhorst-Pack method'©. For all the five 
structures considered in this study, we used unit cells consisting of 8 Si atoms. The con- 
ventional tetragonal cells were used for the diamond and /5Sn structures while for the sh 
and sc structures, the 2x2x2 cells of the conventional primitive cells were used. For the 
hep structure, we took the 2 x 2 x 1 of the conventional primitive cell to make an 8-atom 
unit cell. The Monkhorst-Pack parameters listed in Tab. 1 correspond to those used in the 
respective 1st Brillouin zones of the 8-atom unit cells just described. 

The total energies, i.e. E (apteral , a 2 ), were evaluated at various lattice constants 
^lateral and a z of the unit cells for every 0.05 A. The ranges of lattice constants included 
in the calculations were listed in Tab. 1. The diagonal components of a uniaxial stress, i.e. 
Px = Py an d P z (we shall for convenience refer to as 'pressures'), were obtained from the 
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partial derivatives of the evaluated total energy surface E(ai a terah &z), i-e- 

p _ p _ 1 dE(ai ateral ,a z ) _ 1 dE(ai atera i,a z ) 

X ~ V ~ ^a z da lateral k and A da z |a — (ij 

^lateral 

The A is the area of the lateral unit cell, i.e. af ateral for diamond, /3Sn and sc structures and 
■^afaterai f° r sn an( ^ ^ C P structures. Only structures under pressures of less than 20GPa 
were included in the following discussions. Calculations with a higher energy cutoff and 
different k-point set were considered 0. The change in transition pressure due to numerical 
errors is expected to be less than 0.5 GPa. 



3 Results at hydrostatic condition 

The calculated energy-versus-volume relations for the five structures under hydrostatic 
condition, i.e. P x = P z , were plotted in Fig. 1. The lattice constants of the previously 
defined 8- atom unit cells for the five phases at zero pressure were listed in Tab. 2. They are 
in good agreement with the experimental and previous theoretical results'IH —\. Within 
the pressure ranges we considered in this study, i.e. P x < 20GPa and P z < 20GPa, the 
stable phases were found to be diamond, (3Sn and sh structures. The hydrostatic transition 
pressures (call Phydro hereafter) are 11.4 GPa and 13.7 GPa for the diamond-to-/3Sn and 
/3Sn-to-sh transitions respectively. These are consistent with previous studies^!' iil In 
the present studies, we did not include a phase of space group Imma which have been 
identified, both experimentally'!^' and theoretically'!!!]' HH, to be stable between the /3Sn 
and sh phases. However, the main features of differences in the phase transitions between 
the hydrostatic and uniaxial compression, which are what we would like to address in the 
present work, can already be demonstrated by considering the five phases included here. 
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4 Possible patterns of uniaxial phase transitions 



As discussed in the paper I, the appropriate physical quantity to consider in determining 



the stable phases under uniaxial condition isE-^ 



AH = [(E 2 -E l ) + P av (V 2 - Vi)] + (P z - P x )[(^(Vi ~ V 2 ) + J Uydl z )] (2) 
= AH + AW (3) 

where the limits of the l z in the integral are the lattice constants of the crystal phases 
1 and 2 under uniaxial pressures P x and P z . The internal energies (E) and volumes (V) 
(subscript "1" and "2" for the phases 1 and 2 respectively) are those when the phases are 
under the uniaxial pressures P x and P z while P av is the averaged pressure, i.e. (2P x + P z )/3. 
Phase 2 is relatively more stable than phase 1 if the quantity AH is less than zero and the 
phase-transition pressures are those when the condition AH = is met. In addition to the 
AH term, there is a AW which is not present under hydrostatic condition. The magnitude 
of AW depends on the path of the integration and reduces to zero only under hydrostatic 
condition. The physically most reasonable path for the integral is the shortest path, as 
discussed in paper I, and the integral can be written as 

g(o*, a - a Z}1 )(A 2 + VA 2 Ax + A{) (4) 

for phases considered in the present study. By using equation (4), the quantity AW can 
be rearranged as 

AW =^{P z -P x ){^fT 2 + ^T l ){a z , 2 ^fTi-a z , l ^fT 2 ). (5) 

We shall investigate the possible patterns of uniaxial phase transitions by studying the ah 
initially obtained phase-transition curves of AH = on the P x -P z plane for any two chosen 
phases of the five Si structures, i.e. total ten cases will be studied. The possible patterns 
of AH = will be discussed in terms of the contributions from AH and AW respectively 
and analyzed using their explicit formuli (equations (2) and (5)). 

All the possible patterns of the uniaxial phase-transition curves AH = are 
schematically plotted in Fig. 2a and 2b which correspond to the cases whose Phydro of the 
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two considered phases are positive and negative respectively. For the cases with positive 
Phydro-, if there exist uniaxial phase-transition pressures (both P x and P z ) being smaller 
than the Phydro-, they are said to have lower uniaxial phase-transition pressures. These 
phase-transition curves have positive gradients around Phydro on the P x -P z plane (the thick 
dashed and the thick dotted-dashed lines in Fig. 2a). On the other hand, in the cases of 
higher uniaxial phase-transition pressures, i.e. if all the uniaxial phase-transition pressures 
having either P x or P z or both being larger than the Phydro, the gradients of AH = are 
negative around Phydro on the P x -P z plane (the thick solid line in Fig. 2a). For the cases 
with negative Phydro, if the gradients of AH = are negative around Phydro (the thick 
solid line in Fib. 2b), there is no phase transition between the two considered phases when 
under compression. However, if the gradients of AH = are positive around the negative 
Phydro (the thick dashed and the thick dotted-dashed lines in Fib. 2b), the phase transitions 
between the two considered phases occur only when under uniaxial compression, in spite 
of no phase transition when under hydrostatic compression. 

From Eq.(2), the AH term depends explicitly on the energy and volume differences 
and the averaged pressure. If the energy and volume differences of the two considered phases 
depend little on the pressure anisotropy, i.e. P z — P x , the contours of constant AH are 
expected to be straight lines of negative gradients on the P x -P z plane. This is exactly 
what was found for most of the cases we study here, except for the cases of diamond-sh, 
diamond-sc, /3Sn-hcp and sh-sc whose magnitudes of AH are too tiny to have the constant- 
AH contours of reasonables sizes plotted. Under uniaxial condition, the AH alone can not 
lead to lower phase-transition pressures. 

The other contribution to AH is the AW of equation (5). Ab initio calcula- 
tions showed that the constant-AVF contours on the P x -P z plane are mainly straight lines 
of positive gradients. Whether the contribution of AW is large enough to overbalance 
the negative-gradient constant-AiJ contours and leads to the positive-gradient constant- 
AH contours depends strongly on the magnitude of (a 2j 2V^4i~ — a z,\\T^i) in the AW . We 
found that the zero-pressure structural parameters already give a good estimate of the im- 
portance of AW . The individual values of a z and \fA may change considerably when the 
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phases are under pressures, however, the values of (a z ^y/Ai — a^iV^i) remain roughly a 
constant. In the cases where both considered phases are approximate isotropic structures 
at zero pressure, i.e. \fA ~ a z in Tab. 2, the magnitude of {a z ,2\[M- — a z ,iV^i) W1 U 
be tiny. The data in Tab. 2 suggest that for the phase transitions between diamond and 
sh, diamond and sc, and between sh and sc, the contributions of AW to the AH should 
be small. Our ab initio calculations indeed verify that the gradients of AH = around 
Phydro for these three cases are all negative. In Fig. 3 we present the ab initially obtained 
AH = for the diamond-to-sh (the thick solid line) and diamond-to-sc (the thick dash 
line) phase transitions. These are the two cases having higher uniaxial phase-transition 
pressures. As the Phydro for the phase transition from sc to sh is negative (corresponding 
to the schematically thick solid line in Fig. 2b), the sh phase is always more stable than 
the sc phase under positive external pressures. In general, the sign of Phydro can be decided 
from the zero-pressure properties of the phases as well. For two considered phases, if the 
one which has larger zero-pressure volume but lower zero-pressure energy, the Phydro is 
positive. On the other hand, if the one which has larger zero-presure volume also has larger 
zero-pressure energy, the Phydro is negative. 

In the cases where, at zero pressure, the structure of one phase is very anisotropic 
while that of the other is approximate isotropic, the contribution of AW to AH will be 
significant, due to the component of {a Zj <i\[A\ — a z ,\\fA2) in AW , and lead to positive- 
gradient constant-A^T contours. These include the phase transtions between diamond and 
/3Sn, diamond and hep, /3Sn and sh, /3Sn and sc, /3Sn and hep, sh and hep, and between sc 
and hep. The sign of Phydro results in different patterns of AH = 0. With important con- 
tribution from AW , cases of positive Phydro have lower transition pressures (Fig. 2a) while 
for the cases of negative Phydro the phase transitions occur only at uniaxial condition, i.e. 
Px 7^ Pz (Fig. 2b). In the seven cases considered now, those with positive Phydro values are 
diamond-/3Sn, diamond-hep and /3Sn-sh cases. The ab initially obtained phase-transition 
curves AH = for the three cases (shown in Fig. 4) confirm the conjecture that in these 
three cases, uniaxial compression greatly reduce the phase transition pressures. That is, 
the gradients of AH = for these three cases are all positive around Phydro- Whether 
the AH has gradient larger or smaller than one can be decided by the relative structural 
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parameters of the two considered phases at zero pressure. In the cases of diamond-to- /3Sn 
and diamond-to-hcp transitions, both the /3Sn and hep phases have larger y/A than a z while 
the diamond phase has equal magnitude of y/A and a z at zero pressure. Phase transitions 
from diamond to either (3Sn or hep are expected to be easier if the P z pressure is larger 
than the P x as under this condition the dimensions of the compressed diamond phase are 
forced towards the anisotropic ones (y/A > a z ) as those of /3Sn or hep. Similarly in the 
case of /3Sn-to-sh transition, the transition is expected to be easier if P x is larger than P z 
as the transition is from the anisotropic phase of /3Sn (y/A > a z ) towards the approximate 
isotropic sh phase. These intuitive ideas are supported by the explicit formula of AW in 
Eq. (5). However, this way of thinking should not be pushed too far as the phase transition 
is first order and there is a discontinuity in y/A, a z , as well as the structures. 

The remaining four cases, i.e. /JSn-sc, /3Sn-hcp, sh-hep and sc-hep which have 
important contribution of AW and negative Phydro, should have the AH = contours 
look like either the thick dashed or the thick dotted-dashed line in Fig. 2b. However, 
the Phydro for the /3Sn-hcp is so large in magnitude that the AH = is still out of the 
positive-pressure quardrant of the P x -P z plane we included in our ab initio studies, i.e. in 
the region of P x < 20GPa and P z < 20GPa. For the case of /3Sn and sc phases, as the 
transition from /?Sn to sc is to transit from the phase of y/A > a z to the more isotropic 
phase, the transition will be easier if P x > P z . This is exactly what we obtained from 
ab initio calculations (the thick dotted-dashed line in Fig. 5). For the cases of sh-to-hep 
and sc-to-hep transitions, the situations are opposite to the /3Sn-to-sc transition as now the 
transitions are from the more isotropic phases (sh and sc) to the phase having y/A > a z 
(hep). The phase transition line AH = is expected to have the property of P x < P z . 
In the case of sc-to-hep, the contribution of the AH is tiny and the AW is the dominant 
factor, therefore the AH = phase-transition line is almost parallel to the P x = P z line. 
Note that for the last three considered cases of /3Sn-sc, sh-hep and sc-hep, although there 
is no phase transition at hydrostatic condition, the phases do transit when under uniaxial 
compression, i.e. P x / P z . That is, there might be phase transitions which cannot be 
observed when the phases are under hydrostatic compression but can be observed when 
under uniaxial compression. 
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5 Phase Diagram of Si 



In the previous section, we discuss the possible phase-transition patterns of Si by considering 
only two phases of the five strucutres at a time. However, the phase diagram should be 
constructed when all the phases are taken into consideration. With all the five structures 
taken into consideration, we obtain the phase diagram of Si under uniaxial compression 
as shown in Fig. 6. The stable phases are diamond, /3Sn and sh which are the same as 
those when under hydrostatic condition but stabilize at different regions of the P x -P z plane. 
Under hydrostatic compression, the phase transition is from diamond to (3Sn on increasing 
pressures and to sh if the pressures are increased even further. However, the phase diagram 
suggests that if the external pressures are kept under the condition of P x > P z , on increasing 
pressures, the phase transition should go directly from the diamond to the sh phase (e.g. 
the horizontal dashed arrow in Fig. 6). Similarly, the sh-to-/3Sn transition on increasing 
pressures, which is not possible when the phases are under hydrostatic compression, can be 
observed, according to the phase diagram, if the applied uniaxial pressures are varied from 
the condition of P x > P z to that of P x < P z (e.g. the vertical dashed arrow in Fig.6). 



6 Conclusion 



Based on a previously proposed thermodynamic analysis, we have studied the possible 
phase-transition patterns between phases by carrying out ab initio calculations for five Si 
phases, i.e. diamond, /3Sn, sh, sc and hep. We have shown that it is possible to guess the 
phase-transition pattern of two considered phases by inspecting their structural parameters 
at zero pressure. For cases with positive Phydroi uniaxial compression can greatly reduce 
the phase-transition pressure if the structure of one phase is approximate isotropic while 
the other is very anisotropic. For cases with negative Phydro, it is possible to have phases, 
which do not undergo phase transitions when under hydrostatic compression, do undergo 
phase transitions when under uniaxial compression. With all the five phases taken into 
consideration, the phase diagram of Si when under uniaxial compression was constructed. 
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It was shown that, in spite of the same stable phases were found as those under hydrostatic 
compression, very different phase transitions were possible depending on the conditions of 
the applied uniaxial pressures. 

The author gratefully acknowledges W. H. Huang for invaluable discussions. This 
work wass supported by the National Science Council in Taiwan. The computer resources 
were partly provided by the National Center for High-Performance Computing in HsinChu, 
Taiwan. 
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Tables: 

Table 1. The calculational details used in the present study. See the text in section 2 for 
details. 



diamond /3-Sn sh sc hep 



8-Si unit cell 
k-point sets 

^•lateral (A) 
a z (A) 



tetragonal 

(4 4 4) 
4.60~5.80 
4.60-6.20 



tetragonal 

(8 8 20) 
6.30~6.95 
2.40~2.90 



2x2x2 

(9 9 9) 
4.65~5.40 
4.30~5.35 



2x2x2 
(10 10 10) 
4.50~5.20 
4.50~5.20 



2x2x1 
(7 7 17) 
5.70~7.70 
2.10~3.70 



Table 2. The calculated structural parameters and total energies of the five phases at 
zero pressure. 





diamond 


P-Sn 


sh 


sc 


hep 


^lateral {A) 


5.459 


6.815 


5.295 


5.067 


7.096 


Va (A) 


5.459 


6.815 


4.927 


5.067 


6.603 


a z (1) 


5.459 


2.650 


4.991 


5.067 


3.002 


V (1 3 /Si) 


20.334 


15.385 


15.147 


16.257 


16.363 


E (eV/8Si) 


-43.546 


-40.934 


-40.777 


-40.562 


-40.021 
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Figure Captions: 
Figure 1. 

Energy of the five Si phases, i.e. diamond, (3Sn, sh, sc and hep, plotted versus volume. 
Figure 2. 

Schematic plots showing possible patterns of phase-transition curves (the thick lines) : (a) 
cases with Phydro > and (b) cases with Phydro < 0. The phase-transition curves corre- 
spond to the AH (P x , P z ) = contours. The thin solid lines are the line of P x = P z and the 
thin dashed lines are to guide the eyes to show the hydrostatic phase-transition pressures 

ro- 

Figure 3. 

The ab initially obtained phase-transition curves, i.e. the AH (P X ,P Z ) = contour, for the 
diamond-to-sh (the thick solid line) and diamond-to-sc (the thick dashed line) phase tran- 
sitions. In both cases, the diamond structure stabilizes at the lower and left-handed regions 
of the phase-transition lines. The thin solid line, i.e. the line of P x = P z , corresponds to 
the hydrostatic condition. 
Figure 4. 

The ab initially obtained phase-transition curves for the diamond-to-/3Sn (the thick dashed 
line), diamond-to-hep (the thick dotted-dashed line), and the /3Sn-to-sh (the thick solid 
line) transitions. In both the diamond-/3Sn and diamond-hep cases, the stable phases of 
the lower and right-handed regions of the transition lines are the diamond structure. In the 
/3Sn-sh case, the stable phase of the lower and right-handed region of the transition line is 
the sh structure. The thin solid line corresponds to the hydrostatic condition, i.e. P x = P z . 
Figure 5. 

The ab initially obtained phase-transition curves for the /3Sn-to-sc (the thick dotted-dashed 
line), sh-to-hep (the thick dashed line), and the sc-to-hep (the thick solid line) transitions. 
The stable phases in the lower and right-handed regions of the transition lines are the sc, 
sh and the sc structure for the /3Sn-sc, sh-hep, and sc-hep cases respectively. The thin solid 
line corresponds to the hydrostatic condition, i.e. P x = P z . 
Figure 6. 

The phase diagram of Si under uniaxial compression of P x < 20GPa and P z < 20GPa when 
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the five phases, i.e. diamond, /3Sn, sh, sc, and hep, were taken into consideration. The 
phase boundary between the /3Sn and sh phases is denoted by the dashed line in order to 
show its difference from the hydrostatic line of P x = P z (the thin solid line) . The horizontal 
and vertical arrows show the possible phase transitions which are not possible when under 
hydrostatic compression, i.e. along the thin solid line of P x = P z . 
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